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Abstract 



We construct a solution to the equation of motion of Hamiltonian lattice QCD 
in the strong coupling limit using Wilson fermions which exactly diagonalizes 
the Hamiltonian to second order in the field operators. This solution obeys 
the free lattice Dirac equation with a dynamical mass which is identified with 
the gap. The equation determining this gap is derived and it is found that 
the dynamical quark mass is a constant to lowest order in N c but becomes 
momentum dependent once 0{\/N c ) corrections are taken into account. We 
interpret our solution within the framework of the N-quantum approach to 
quantum field theory and discuss how our formalism may be systematically 
extended to study bound states at finite temperature and chemical potential. 
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I. INTRODUCTION 



The strong coupling approximation has played an important role in the development of 
QCD lattice gauge theory from its very inception. In the renowned paper by Wilson this 
approximation was invoked to demonstrate quark confinement on the Euclidean space-time 
lattice. Soon thereafter Kogut and Susskind formulated the Hamiltonian lattice gauge 
theory and concluded that in the strong coupling limit the quark dynamics is best described 
by a collection of non-Abelian electric flux tubes with quarks attached at their ends. This 
was followed by the work of Baluni and Willemsen who used a variant of the Kogut-Susskind 
formalism to demonstrate quantitatively that dynamical chiral symmetry breaking indeed 
takes places in lattice QCD at strong coupling ||. Finally, calculations by Kogut, Pearson 
and Shigemitsu || and by Creutz || suggesting the absence of a phase transition between 
the strong and weak coupling regimes of QCD motivated numerous studies in this subject 
which continues to date. 

The Hamiltonian formulation of strong coupling QCD using Wilson fermions was first 
studied in detail by Smit || who derived an effective Hamiltonian using the 1/N C approx- 
imation and applied the spin wave theory of magnetism to determine the vacuum energy 
density and the excitation spectrum. In the absence of the current quark mass and the 
Wilson term the theory posseses a U (4/Vf) symmetry. Smit has shown how this symmetry 
is spontaneously broken down to U(2N{) <8> U(2N{) accompanied by the appearance of 8iV f 2 
Goldstone bosons. A finite current quark mass also breaks the original U(4Nf) symmetry 
explicitly down to U(2N{) <g) U(2N{). Introduction of the Wilson term explicitly breaks the 
latter symmetry further down to U(Nf) solving the fermion doubling problem. In addition 
Smit has proposed how quantum corrections appearing at 0(1/N C ) may be used to resolve 
the Ua(1) problem. Le Yaouanc, Oliver, Pene and Raynal 0] later reexamined Smit's work 
by applying the Bogoliubov-Valatin variational method to the same effective Hamiltonian 
and found that the vacuum is not chirally degenerate and concluded that the massless pseu- 
doscalar boson found by Smit is not a Goldstone boson. 

In this letter we report on a new approach to study strongly coupled Hamiltonian lattice 
QCD using Wilson fermions. Our approach differs from those adopted in |5| and [7| in 
that we explictly construct a solution to the field equations of motion derived from the 
Hamiltonian. This solution exactly diagonalizes the Hamiltonian to second order in field 
operators and obeys the free lattice Dirac equation where the mass plays the role of the gap. 
Using the equation of motion we derive the gap equation for the dynamical quark mass to 
0(1/N C ). Our results for the dynamical quark mass is qualitatively different from the one 
found in J7|. 

With our solution the interpretation of the chiral condensate being the order parameter 
becomes transparent as well as the observation that the broken symmetry phase is the 
energetically favored phase. In addition, our approach admits to a systematic extension 
to the study of bound states both in free space and at finite temperature and chemical 
potential. This is accomplished by interpreting our solution within the context of the N- 
quantum approach to quantum field theory |8||| which we shall discuss in the concluding 
section. 
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II. NON INTERACTING WILSON FERMIONS 



We begin with a brief review of the free lattice Dirac Hamiltonian with the Wilson term. 
Although elementary, we shall see that the results presented here will play crucial roles in 
the construction of a solution of the strongly coupled QCD. We adopt the notation of Smit 
H with unit lattice spacing (a = 1) and temporarily suppress color and flavor indices. Then 
the free lattice Dirac Hamiltonian is 



H° = ^.Y. [^(x)j ji^(x + ni) -^ f (f +n ; )7o7^(^)] +M^^ t (f) 7o ^ 

— x,l 



[X 



- 1 x,l 



(2.1) 



where the third term is the Wilson term with < r < 1. For r = there is an eightfold 
fermion multiplicity which is removed when r ^0. 

At each lattice site the free Dirac field in configuration space is given by 



(2.2) 



with fi denoting the Dirac index. The excitation energy u(p) will be determined shortly. 
The annihilation operators b and d annihilate the non-interacting vacuum state |0). For 
our purpose it is not necessary to know the structure of the spinors £ and 77. The only 
assumption that we shall make is that the creation and annihilation operators obey the free 
fermion anti-commutation relations 



tf(j>),b{q)] t = \d\p),d(q) 



(2.3) 



Using this assumption we can recover the anti-commutation relations for the field operators 



provided that £ and 77 satisfy the relation 



We normalize the spinors by demanding that the number density is given by 
M = ]T : ¥(t, x )tf (t, x ) : = 2 £ (b\p )b{p ) - d\p )d{p)) 



(2.4) 



(2.5) 



(2.6) 



where the symbol : : denotes normal ordering and the factor of 2 accounts for the spin 
degrees of freedom. Eq. ( |2.6| ) fixes the normalizations of £ and 77 to be 

UP)&(P) = %(pH(p) = 2 (2.7) 

= vI(p)U-p) = (2-8) 

which are consistent with Eq. 
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We now go over into momentum space where we perform all our calculations. In mo- 
mentum space the charge conjugaton symmetric form of H° is [] 



p I 

where the momentum dependent mass term is given by 

M{p) = M — r V] cos(p • hi) 
i 

The free Dirac field becomes 

which is used to derive the equation of motion corresponding to Eq. 
M{t,p) = : [*(t,p),H°]_: 

^sinQo- hi)j li + M(p)7o^(t,p) 



1/2 



From Eq. ( gl3|) one obtains the excitation energy 

^(p) = (E sin2 (p • "0 + m2 (p) 

and the equations of motion for the £ and 77 spinors 

w(p) £(p) = (^sin(p • n«)7o7i + M (^70) £(p ) 
w(p) r?(-p) = -(^sinQo- r^)7o7« + M(p)j Jr)(-p] 



(2.9) 
(2.10) 

(2.11) 

(2.12) 
(2.13) 



(2.14) 



(2.15) 
(2.16) 



When r = these equations of motion are relativistic near the eight corners of the 
Brillouin zone denoted by n = (0,0,0), 7r x = (tt, 0,0), n y = (0, 7r, 0), 7r z = (0,0, tt), 7r xy = 
(7r, 7r, 0), 7f xz = (7r, 0, 7r), 7r yz = (0, 7T, 7r) and 7T xyz = (ff, 7T, 7r). The excitation energies at these 
values of momenta are equal which corresponds to the eightfold multiplicity mentioned 
above. This degeneracy is lifted when r ^ due to the momentum dependent mass term 
Eq. ( |2.10| ). Using the equations of motion for £ and 77 it is a simple excercise to show that 
the off-diagonal Hamiltonian vanishes and that the vacuum energy is given by (0|ff°|0) = 
— 2V Yjpoj(p). Finally, we construct positive and negative energy projection operators A + (p ) 
and A~(p) as follows 



x Our convention for the Fourier transform from configuration to momentum space is 9(x 
J2fi^(p)e p ' x which implies that the volume V is given by V = J2g = 8pg. 
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a+(p ) = e(p ) ® ^(p ) 
i 

~ 2 



1 \ M(p) 

w(p)Y w (p) 



A (p) = r/(— p) <S> ry{— p) 
1 



2^ sin(p ■ n,)7o7/ t^7o 



(2.17) 



(2.18) 



These projection operators will be used extensively below. Note that they obey the condition 



A+(p) + A - (p) 



a/3 



5, 



a/3 



(2.19) 



as is required by Eq. 



III. EFFECTIVE HAMILTON! AN AND THE ANSATZ 

For the sake of comparison we use the same effective Hamiltonian obtained by Smit in 
Eq. (3.13) of || which is also used in [JJ]]. The charge conjugation symmetric form of Smit's 
Hamiltonian in momentum space is 



Pl,P2 



K 



( Vi+-+p 4 ,o 

pi,...,P4 £ 



gi((^i+^)-ni) _|_ e i((p3+P4)-"0 



(Sl)^ (^) 7 (ft )(^),(P4) - (S/) 7 , (* W ),(ft)(*i,) 7 (p4) 



(3.1) 



where S ; = —i^oli — * r 7o)- We denote color, flavor and Dirac indices by (ab), (a/3) and 
(fjjWyS), respectively. Summation convention for repeated indices is implied. The three 
parameters in this theory are the Wilson parameter r, the current quark mass Mq and 
the effective coupling constant K which behaves as 1/g 2 with g being the QCD coupling 
constant. 

Our ansatz which diagonalizes Eq. (|3.1| ) to second order in the field operators is decep- 
tively simple. It is the solution of the free lattice Dirac equation Eq. (|2.13| ) with r = and is 
given by Eq. ( |2.11| ) appropriately modified to include color and flavor degrees of freedom.^ 
However, the mass term M(p ) in Eq. (|2T3| ) is now the dynamical and not the current quark 
mass. We shall show below that this dynamical mass is the solution of the gap equation. 
In addition, the annihilation operators b and d in our ansatz now annihilate an interacting 
vacuum state \Q) and not |0). Since we shall be working only in the space of quantum 



2 This modification will of course affect the anti-commutation relation Eq. (|2.3|). 
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field operators there is no need to specify the structure of \Q). Aside from these non-trivial 
modifications our ansatz obeys all the properties of the free Dirac field described above. For 
example, with this ansatz the chiral condensate is given by 



TV Nt 

{G\m<3) = "V E [ Tr (7oA-(p)) " Tr( 7o A + (p)) 



-2N c N f J2 



Mjp) 
uj(p) 



(3.2) 



Hence the chiral condensate is determined by the dynamical quark mass and clearly plays 
the role of the order parameter. The basic idea of our approach is to use the equation of 
motion for H c r to derive and solve the gap equation for the dynamical mass. 

Using our ansatz we proceed to derive the vacuum energy density and the second order 
off-diagonal Hamiltonian by renormal ordering H c r with respect to \Q). The former is given 
by terms involving two contractions while the latter is derived by retaining terms involving 
a single contraction. Our results for the vacuum energy density is 

^(Q\H cS \g) = -2M 7V c ^Tr [A + (p) 7( / 

p 

-^ee a ;(p)^(?) 

p,q I 



<S> |cos (p — q) - hi 

+ cos {p + q ) • hi 



( S /) jUI /( S i) 7 5 + ( S ')j^( S i) 7 <5 
S Mi ,( S 7 <5 + ( S 0^( S 7 5 



+ cos (p + q ) - hi 
while the second order off-diagonal Hamiltonian is found to be 
H oS \G) = |MoE(7o) M 4(g)^(-g) 

+/eea;(?)$?)^(-?) 



cos {p - q ) ■ h x (W^ls + (E,)^(S,)^ 



+ cos (p + q ) • n, ( (^i)l(^)l s + (£*U^) 75 



EE [2A+ (?) - ^] C(?M-<?) 



JV C 4 - - , 

p,q I 



(3.3) 



6 



+ cos (p + q)-h t {{^iMis + (^W^) 7 s) }« 



®bUqKa(-Q)\G) 



(3.4) 



We see from Eq. Q3.4p that the excitation spectrum of the effective Hamiltonian involves 
color singlet quark-anti-quark mesonic excitations coupled to zero total three momentum. 



IV. EQUATION OF MOTION AND DIAGONALIZATION OF H EFF 



We shall now demonstrate that H q vanishes exactly by exploiting the equation of motion 
for H e s . For this purpose it suffices to consider the equation of motion for the up quark field 
u(t, q) given by 



K^a) u {t,q) = ■ (u a ) u (t,q),H ( 



'cff 



M (70 ; 



/j,S 



(4.1) 
(4.2) 



c p 1 



cos (p- q)-nt [(^l\ v ^l)U + (SO^CSi)^) 
+ cos (p + q)-n t [fr^Wh + (Si)^(Si; 



1 K 



EE[2A+(?)-^ 



iV c 4 . _ , 

P>9 ' 



+ cos (p + q) ■ nil ^i)U^i)\s + (E)^(E) 7 , 



(u a ) s (t,q) (4.3) 



The crucial step in diagonalizing H e g is to recall that u a is also the solution of the free lattice 
Dirac equation with r = and therefore one can identify Eq. Q4.1|) with Eq. ( |2.13| ). 



u a , H eS 



a ai 11 r =0 



(4.4) 



Then the second order off-diagonal Hamitonian vanishes due to the equation of motion for 
r\ Eq. (|2.16|) and the orthogonality condition Eq. fl2.8|). 



Hoz\Q) = ?(q) f E sin(g • n,)7o7* + M(q ) 7 oj ^(-^c^^^^C-^)^) 

= (4.5) 
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We can also use Eq. ( |4.4|) to simplify the vacuum energy density to the following form. 

^(G\H e n\Q) = -N c £ {3^(1 + r 2 ) + 2u(q) + j (4.6) 

Therefore, the difference in the vacuum energy densities between the phase with M(q ) = 
and M(q) 7^ is given by 

= —(G\H cS \Q)\ M (q)=o - y(G\H cS \g)\ M (q)^o 



^M + 2^sin 2 (g-n i ) + M 2 (g)^ ' -2^sin 2 (g-^)) ' 



> (4.7) 

Hence the phase with a finite gap is the energetically favored phase. It remains to derive 
the gap equation and to solve it for the dynamical quark mass. 



V. THE GAP EQUATION AND THE DYNAMICAL QUARK MASS 



The gap equation is obtained by using Eq. ( |4.4j ) and equating the coefficients of the 70 
operator. Using Eq. ( p.l?| ) we find 

M(^M„ + ^(l- r *)£^M 

1 M(p) f 2 1 2 1 

Thus the dynamical quark mass is determined by a three dimensional non-linear integral 
equation and is momentum dependent. This dependence comes from the 0(1/N C ) correction 
to the gap equation. 

The solution to the gap equation to lowest order in iV c is shown in Figure la for the 
case of M = and r = 0. We see that the dynamical symmetry breaking takes place only 
above a critical coupling constant of about Kq 0.73. From Eq. ( |5.1|) we see that K c must 
increase as r is increased in order for the gap equation to have a solution and this situation 
is shown in Figure lb. In fact, the critical coupling constant approaches infinity as r — > 1 
and the assumption of strong coupling (small K) breaks down in this limit. 

However, when r = 1 there is no 0(N®) contribution to Eq. (|5.1|) and the dynamical 
quark mass may be written as 

M(q)=M + B(K)J2cos(q-h l ) (5.2) 

assuming that M(—q) = M(q). In this case the solution to the gap equation is shown in 
Figure 2 again for the M = case. At this value of the Wilson parameter we find solutions 
to the gap equation only above a value of Kq w 0.84. However, the dynamical quark mass 
is one order of magnitude less than that obtained with only the 0(N®) contribution. This 
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result suggests that the momentum dependence of the dynamical quark mass may be treated 
correction to the dominant constant term. 
Our results for the dynamical mass is qualitatively different to the one obtained by 
Le Yaouanc, Oliver, Pene and Raynal in Eq. (6.6) of 0. Not only their dynamical mass 
is momentum independent, but the chiral symmetry breaking takes place for all values of 
the coupling constant K > 0. In addition, from Eq. (|3.2|) we see that in the Mq — > limit 
the chiral condensate is finite only above a certain value of the critical coupling constant 
Kq for the examples presented here. Above Kq, the chiral condensate is not a constant 
but a function of K. This situation is in contrast to the behaviour of the chiral condensate 
obtained in |7[ where it is always a finite constant in the broken symmetry mode. 



VI. CONCLUSION AND OUTLOOK 

In this work we have constructed a simple solution to Smit's effective Hamiltonian for the 
strongly coupled lattice QCD. Our solution obeys the free lattice Dirac equation with a dy- 
namical quark mass. This mass is the solution of the gap equation which is determined from 
the equation of motion. The creation and annihilation operators for quarks and anti-quarks 
obey the well-defined free fermion anti-commuation relations and annihilate an interacting 
vacuum state. Since quantites of interest are solely determined by the anti-commutation 
relations there is no need to specify the structure of this state. 

The solution presented here exactly diagonalizes the Hamiltonian to second order in 
the quark field operators. We find that the elementary excitations of the theory consist of 
quark-anti-quark color singlet states coupled to zero total three momentum and that the 
phase with broken chiral symmetry is the energetically favored phase. Morever in this phase 
the gap of the theory is found to depend both on momentum and on the effective coupling 
constant. 

However, our solution is by no means complete. It may be systematically improved 
to calculate the masses, widths and coupling constants of all the bound states allowed by 
the Hamiltonian if our approach is to be interpreted within the context of the N-quantum 
approach (NQA) to quantum field theory MM as follows. NQA is a method to solve the 
field equations of motion by expanding the Heisenberg fields in terms of asymptotic fields 
obeying the free field equations of motion. This expansion, which varies from theory to 



theory, is known as the Haag expansion 1 10 ] 



The first term in the Haag expansion is simply a single free asymptotic field. Previous 
applications of the NQA to various field theories have shown that the use of this first 
order term alone can reproduce the known mean field results not only in free space but 



also at finite temperature (T) and chemical potential (/i) H|1T|. These additional variables 
are introduced by subjecting each of the field operators in the Haag expansion to thermal 
Bogoliubov transformations just as in thermal field dynamics. The solution presented in 
this work can be interpreted as the first order Haag expansion of the Heisenberg quark 
fields in the effective strong coupling QCD Hamiltonian. Since the time coordinate is not 
discreticized in the Hamiltonian (Kogut-Susskind) formulation of lattice gauge theory one 
can naturally introduce and work with the concept of an asymptotic field. 

We can improve our solution by taking the second order terms in the Haag expansion 
into account. In our case each of these terms would be a product of a fermionic (quark) and 
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a bosonic field, the latter corresponding to the various bound states allowed by the Hamil- 
tonian. They are the elementary excitations identified in this work. The coefficients of the 
second order terms are interpreted as amplitudes for formations of these bound states. Self- 
consistency equations for these bound state amplitudes are derived by using the equation 
of motion. These bound state amplitudes have the same number of kinematical variables as 
in non-relativistic theories and are independent of Bethe-Salpter amplitudes. The idea of 
the NQA is to solve for these amplitudes and thus construct a solution to the field theory 
beyond the mean field approximation. 

The program described above have recently been carried out for the effective instanton 
induced 't Hooft interaction at finite T and \i In constructing a solution to the 't Hooft 
model beyond the mean field theory the masses, widths and coupling constants of the a 
and diquark states have been determined. It would be interesting to apply the second order 
NQA to the Hamiltonian considered in this work and to calculate the properties of bound 
states of the theory and to construct an improved solution to the strong coupling QCD 
Hamiltonian at any T and \x. This solution can then be used to obtain a qualitative picture 
of the QCD equation of state. 
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FIGURE CAPTIONS 



FIGURE 1. The solution to the gap equation Eq. ( |5.1| ) to lowest order in N c for M = 0. 
(a) Dynamical quark mass m for r = as a function of the effective coupling constant K. 
The critical coupling constant Kq is approximately 0.73. (b) Kc as a function of the Wilson 
parameter r. 

FIGURE 2. The solution to the gap equation Eq. ( |5.1| ) f° r r = 1 and M = 0. The 
coefficient of the momentum dependent dynamical quark mass B from Eq. (|5.2p is shown 
as a function of the effective coupling constant K. The critical coupling constant Kc is 
approximately 0.84. 
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